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CN . Abstract 

Q , In antiferromagnetically coupled superlattices grown on (001) faces of cubic substrates, e.g. based 
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on materials combinations as Co/Cu, Fe/Si, Co/Cr, or Fe/Cr, the magnetic states evolve under 
competing influence of bilinear and biquadratic exchange interactions, surface-enhanced four-fold 
in-plane anisotropy, and specific finite-size effects. Using phenomenological (micromagnetic) the- 
' ory, a comprehensive survey of the magnetic states and reorientation transitions has been carried 

|: 

. plane. In two-layer systems [N = 2) the phase diagrams in dependence on components of the ap- 
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plied field in the plane include "swallow-tail" type regions of (metastable) multistate co-existence 
and a number of continuous and discontinuous reorientation transitions induced by radial and 



^ ■ transversal components of the applied field. In multilayers {N > 4) noncollinear states are spa- 

o '. 

Q , tially inhomogeneous with magnetization varying across the multilayer stack. For weak four-fold 

anisotropy the magnetic states under influence of an applied field evolve by a complex continuous 

>< ' 

■ reorientation into the saturated state. At higher anisotropy they transform into various inhomo- 

Cd . 

geneous and asymmetric structures. The discontinuous transitions between the magnetic states in 
these two-layers and multilayers are characterized by broad ranges of multi-phase coexistence of 
the (metastable) states and give rise to specific transitional domain structures. 
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I. INTRODUCTION 



Multilayers built from ferromagnetic layers with various spacer layers include a wide va- 
riety of magnetic film systems that have been intensively studied during last years. Due 
to remarkable phenomena as giant magnetoresistance, exchange-spring behaviour and/or 
exchange bias, and surface enhanced magnetic anisotropy, such nanostructures have already 
found a number of applications and are considered as promising candidates for nonvolatile 
magnetic recording media. ^ On the other hand, nanoscale superlattices and similar struc- 
tures provide convenient model systems to study different aspects of surface magnetism and 
magnetic ordering in confining geometry. 

In particular, much attention has been ^ive^n_to_jnultilayers composed of antiferromagnet- 



ically coupled ferromagnetic nanolayers. 
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Such layered synthetic antiferromagnets can be separated into two classes according to the 
symmetries ruling their magnetic properties: Superlattices with relatively strong uniaxial 
anisotropics include low-symmetry multilayers with effective uniaxial magnetic anisotropy 
in the layer planes, e.g. epitaxial systems deposited on (110), (211) faces of cubic 
^ Also multilayer systems with perpendicular anisotropy belong to this 



substrates 
class. 



m 




Magnetization processes in these nanostructures are strongly influenced by the 
uniaxial anisotropy which is responsible for specific phenomena such as "surface spin-flop" 
^,0] or field-induced metamagnetic jumps. [5] On the other hand, the uniaxial anisotropy 
may be absent in layered systems with higher symmetry. These represent another large 
and intensively investigated class of synthetic antiferromagnetic nanostructures. Superlat- 

cubic substrates, e.g. multilayers 



tices with planar magnetization grown on (001) faces o 
rom materials combinations as Co/Cu l9|, Fe/Si ] 
151 Il6l . Il7| belong to this class. In the case of weak four-fold anisotropy, their magnetic 
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IJ or Fe/Cr, 



properties are mostly determined by the interlayer exchange interactions which may in- 
clude an important biquadratic contribution. ^,11, 3, Evidence of strong biquadratic 
exchange interaction has been given in a number of experimental papers for Fe/Cr two- 
layers and multilayers 13,13, fo'^ Fe/Si(001) multilayers |lll)- Strong effective four-fold 
anisotropics have been found in systems such as Co/Cr(001) or Fe/Cr(001). P. Iisl. Il^ lisj] 
The competition between magnetic anisotropy, applied fields and exchange energies may 
cause complicated magnetic effects and processes. In fact, a great number of novel mag- 
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netic configurations and remarkable reorientational effects in external fields have been found 



in such superlattices. jol IG, 15, 17 1 In particular, recent experimental results using mod- 
ern depth-resolving techniques reveal spatially inhomogeneous magnetic structures, e.g. in 
Fe/Cr(001) superlattice s. |l7l| and specific reorientation effects imposed by four-fold planar 
(tetragonal) anisotropy. Il5| The understanding and interpretation of the complex magneti- 



zation processes found in such systems 
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|20| requires a theoretical underpinning. 
Theoretical activity in this field is largely based on analytical and numerical calcula- 
tio. .0.1. witHin pHe_,o,ca. approach.. ^2^^^ TW .ud.e. Kave 

demonstrated the general validity of the phenomenological models to describe the magne- 
tization processes in antiferromagnetic nanostructures. 

mm 

For the system under 

discussion the phenomenological theory has been developed to describe effects of 90 de gree 
couplings and a concomitant complex evolution of domain structures in Fe/Cr/Fe layers. 
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Four-fold anisotropy effects have been theoretically investigated in Ref. 
structures with N=2 ferromagnetic layers coupled through a spacer. Multilayer systems 
provide also experimental models to study effects of the confining surfaces on antiferromag- 
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netic structures, 
theoretical models in Ref. 



lis context, specific inhomogeneous magnetic states described for 



231 have recently been observed in Fe/Cr superlattices. jl7| The 



existing theoretical results, however, are restricted to special cases and are not sufficient 
for an exhaustive description of the rnagnetic states and field-induced reorientional effects 
observed in recent experiment. [l^Ill, 12, 1^ 14, 1^ 3, 1^ 

In this paper we provide a theoretical analysis of magnetic states and magnetization 
processes in planar antiferromagnetic superlattices with and without four-fold anisotropy 
in magnetic fields applied within the multilayer plane. In such magnetic superlattices the 
exchange interlayer coupling is an oscillatory function of the spacer thickness. Q, Q De- 
pending on the spacer thickness an alternating sequence of ferromagnetic and antiferromag- 
netic interlayer couplings is realized, and, by adjusting the spacer thicknesses, very different 
strengths of antiferromagnetic coupling can be realized. On the other hand, the four-fold 
anisotropy includes interface-induced contributions which implies a strong dependency of 
the effective average anisotropy of each ferromagnetic layer on the layer thickness. [3] Thus, 
the effective magnetic interactions can vary in extremely broad ranges for such multilayers 
in dependence on the chosen materials combinations and thicknesses, see e.g. Ref. where 
for a Co/Cu(001) two-layer system with a wedged spacer layer the ratio between four-fold 
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anisotropy energy and the exchange couphng is changed by orders of magnitude. In con- 
trast to bulk planar antiferromagnets, where an essentially fixed hierarchy for the strengths 
of the magnetic interactions holds, [29[ in these artificial antiferromagnetic systems the ra- 
tios between different magnetic energies, respectively the phenomenological parameters in 
the magnetic free-energy, may assume practically arbitrary values. Moreover, as the in- 
terlayer exchange is weak compared to direct exchange interactions, the fields to induce 
spin-reorientation phenomena are similarly weak and experimentally accessible. 

The rich phase diagrams for these systems precludes an analysis in all details. The 
phase space in terms of the phenomenological parameters includes a large variety of dif- 
ferent magnetic states with a corresponding multitude of spontaneous and field-induced 
phase-transformations. In a first step to such an analysis, all laterally homogeneous states 
in such multilayers must be found. They are the building blocks for a domain theory. 29 



For the case of laterally homogeneous states of each ferromagnetic sublayer one has a sys- 
tem that behaves like an antiferromagnetically coupled chain of Stoner-Wohlfarth particles. 
This simplified one- dimensional model for the behaviour across the multilayer stack also 
yields the limiting case for the magnetization processes with maximum hystereses. Again, 
a direct analysis of all magnetic states even for these one-dimensional models yields an in- 
tricate succession of phase diagrams and magnetization curves. In this paper, we avoid 
the cumbersome task of listing and classifying all solutions and transitions. Instead, we 
provide a broad physically intuitive picture of the physical effects due to the different ex- 
change or anisotropic forces, and those imposed by the confining geometry of the system. 
To this end, we study limiting cases of the model. This includes the case of strictly zero 
anisotropy with and without biquadratic exchange, and the case of infinite anisotropy with 
fixed four-fold orientation of the magnetizations in each layer. For the antiferromagnetic 
two-layer systems (A^=2), i.e. the sandwich structures (experimentally realized as ferromag- 
netic/spacer/ferromagnetic trilayers), we provide a detailed investigation of the magnetic 
phase-diagram for arbitrary orientation of fields in the layer planes. Based on this, we can 
understand the basic magnetic configurations in the multilayers, and we can give a map of 
the topologically different types of magnetic phase diagrams. 

We use standard methods to analyse magnetic phases and transitions within the phe- 



nomenological approach and the theory of phase transitions. 



3l| The one-dimensional chain 



models are considered as composite order-parameters with many components (A^ compo- 
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nents in a multilayer stack composed of ferromagnetic/spacer bilayers) and a character- 
istic structure of couplings between the components defined by interlayer exchange and the 
surfaces. From this point of view, the very rich phase-diagrams and correspondingly com- 
plex sequences of magnetic configurations can be understood. For the general cases of the 
model, the equations for equilibrium and phase stability can be solved only by numerical 
methods. With the methods and results expounded below, one can extend the analysis to 
specific experimental cases in all detail. 

The solutions for the one-dimensional chain models include various field-induced canted 
and inhomogeneous states with a net magnetization. Based on the phase diagrams of these 
models, the evolution of laterally inhomogeneous (domain) states and magnetization pro- 

In this connection, the coexistence regions of different phases 



cesses can be discussed. 
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in the vicinity of discontinuous (first-order) magnetic phase transitions are important. In ex- 
ternal fields, thermodynamically stable domain-configuration from these competing phases 
can be established in extended multilayers This is crucial for an understanding of the hys- 
teretic magnetization processes under coercivity mechanisms. 

The paper is organized as follows. After introducing the model and mathematical tools 
(Sec. II) we consider the effects of the bilinear and biquadratic exchange interactions in the 
next section (Sec. III). In Sec. IV we investigate in detail four-fold anisotropy effects in 
antiferromagnetically coupled two-layers and then discuss the generalization of these finding 
to the case of multilayers. In Sec. IV we discuss domain states and magnetization pro- 
cesses by using quahtative arguments. In the concluding part we discuss possible extensions 
of the theory, and we suggest some useful experiments to enhance our understanding of 
magnetization processes in antiferromagnetic superlattices. 



II. THE MICROMAGNETIC ENERGY 



Let us consider a stack of ferromagnetic plates infinite in x- and y-directions and with 
finite thickness along 2;-axis. The magnetization of the layers is Mj, and there are indirect 
interlayer-exchange couplings through spaces between them. The phenomenological energy 
of the system can be written in the following form 
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= ^ ^Ji nii ■ rrii+i + Jj (nii • rrii+i)^ - H ■ ^ duii - i^^Ki (m^ + mj) (1) 

i=l 1=1 i=l 

where iiij = Mj/Mq*^ (^o*^ = |Mj|) are unity vectors along the i-th layer magnetization. 
Q = Mq '' /Mq designate deviations of the magnetization in the i-th layer from the averaged 
value Mq. We assume that miz=0 , i.e. the layer magnetizations are restricted to the layer 
plane. Jj and Jj are constants of bilinear and biquadratic exchange interactions, respec- 
tively; Ki are constants of the in-plane four-fold anisotropy. The functiona l (ITll ge neralizes 
similar models considered earlier in a number of studies on exchange 16^1211. and 



anisotropy [9 



15 



22 



29| effects in planar antiferromagnetic systems. Within this approach 



the ferromagnetic layers are considered as homogeneously magnetized blocks with constant 
values of the magnetic interactions. This assumption deserves some comment. It is well- 
established that in magnetic nanostructures surface/interface exchange and relativistic inter- 
actions strongly modify electronic and magnetic properties within all volume of the magnetic 
constituents [33]. This means that the values of the exchange or anisotropy parameters, and 
the magnetizations include large interface/surface-induced components which may strongly 
vary across the thickness of the individual layers. However, the hypothesis of magnetic 

homogeneity in the ferromagnetic nanolayers in the models of type has a solid physical 
basis and is justified by successful applications of these models to describe magnetization 
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processes in layered ferro- and antiferromagnetic nanostructures. 
This relies on the fact that in ferromagnetic nanolayers the intrinsic (direct) exchange cou- 
pling are usually very strong and overcome surface/interface induced interactions. Thus, 
they play the dominating role to determine the magnetic order within the layers, which 
reacts also very stiffly on all external and induced magnetic forces. Furthermore, in these 
planar systems the stray field forces confine the magnetization of the layers into their plane. 
As a result, in most of these systems the ferromagnetic layers preserve essentially homo- 
geneous in-plane magnetized states even under the influence of strongly inhomogeneous 
induced interactions. 0. [3^ In this connection it is important to stress here that the ma- 
terials parameters in the phenomenological model should be treated as averages over a 
multilayer period. They comprise (in integral form) all intrinsic and induced energy contri- 
butions of the magnetic states in the ferromagnetic layers. In particular, for the systems 
under consideration they may also include a contribution from magnetism of the spacer 



layers. Thus, in contrast to their bulk counterparts which are considered as constants of 
the magnetic material, in nanoscale multilayer systems the phenomenological parameters 
strongly depend on many physical and geometrical factors and may considerably vary from 
sample to sample. 

Parity of also plays an important role. The model (P) with even number of ferromag- 
netic layers represents systems with fully compensated magnetization. Such superlattices 
can be treated as analogues to bulk collinear antiferromagnets. Superlattices with odd num- 
bers of layers or with unequal thicknesses of layers own a noncompensated magnetization 
which strongly determines their global magnetic properties. In their response to an applied 
field, these structures are similar to bulk ferrimagnets. They could be studied by similar 
methods as used below, but have to be considered as separate class of systems. Concentrat- 
ing on the properties of systems with fully compensated total magnetization, we consider 
only superlattices with even N and equal magnetization in all layers (C? = !)• In this con- 
tribution dedicated to general properties of these antiferromagnetic superlattices, we avoid 
secondary effects which are related to their micro-structure and interface properties, such 
as strains, chemical intermixing etc. With that in view we will study the model (0) mainly 
for the case of identical ferromagnetic layers assuming that the exchange and anisotropy 
constants are equal for all the layers, Ji = J, Ji = J, and Ki = K. Some of the analytical 
results can be generalized to the isotropic (i. e. with = 0) system with mirror symmetry 
about the center of the layer stack, i.e. Ji = Jn-i and Ji = Jn-i- 

In our problem the magnetization of the i'th layer is confined in planar orientations and 
can be described by the angle 9i between the vector mj and x-axis. Thus, calculations of 
the magnetic states for the model (0) can be reduced to the optimization of the function 
$Ar(^i,6'2, ...On)- Following the theory of bulk antiferromagnetism it is convenient to intro- 
duce here the vectors of the total magnetization m and the staggered (or antiferromagnetism) 
vector 1 

N N 

m = J]m,,, l = ^(-l)(^+iVi. (2) 

1=1 i=l 

The energy (PJ is invariant under transformation 1 —1, and, therefore, all magnetic states 
in this model are degenerate with respect to the sign of the staggered vectors. In the following 
only solutions with a definite sign of the staggered vector will be discussed. However, one 
should keep in mind that the full set of solutions includes also those with opposite sign of 1. 
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The magnetic states with opposite sign of 1 behave identically in an applied magnetic field. 
Thus, for the magnetization processes there is no need to distinguish these solutions. 

For the particular case of a two-layer {N = 2) the energy ([Q) can be transformed to the 
following form 

K 3K 
$2 = Ji COS 20 + Ji cos^ 20 - 2H cos cos(6' - i/j) - — cos 40 cos AO — (3) 

where 9 = {9i + 62) /2, (p = {61 — 62) /I and ip is the angle between the H and the x-axis. 

This case of a two-layer is of special interest. The energy (jHl) functionally coincides with 
that of a bulk two-sublattice antiferromagnet. This offers the opportunity for useful physical 
relations and analogies with bulk antiferromagnet ism. On the other hand, the two-layers 
represent the simplest model of antiferromagnetically coupled layered nanostructures. 

The equations minimizing the function $7v(^i; ^2? ■■■(^n) have strongly nonlinear character 
and no analytical solution are available generally. The body of our results have been obtained 
by numerical methods. We solve the coupled equations for equilibria {dW/dOi = 0}i=i,,,N 
by an efficient conjugate gradient minimization. 36| For systems with large anisotropy the 
configurations with the magnetizations oriented in the four-fold easy-axis directions are close 
to mostly metastable minima of the energy. These configurations have been systematically 
searched for global absolute stability and checked for stability. Numerically, this is feasible 
for not too large with present day computers. Stability limits and phase-transitions are 
determined from the evolution of the smallest eigenvalue eo{H,K) of the stability matrix 
{d'^W/d9id9j),i, j = 1 . . . N under changing parameters in ((T)), i.e. ratios J/ J, K/ J, or the 
applied magnetic field. 



III. MAGNETIC STRUCTURES IN ISOTROPIC MULTILAYERS. SPATIALLY 
INHOMOGENEOUS SPIN-FLOP STATES 

A. Exchange interactions in layered antiferromagnets 

We start the investigations of the static configurations minimizing energy (0) from the 
isotropic case (i.e. with Ki = 0). Depending on relations between bilinear and biquadratic 
constants different coUinear and noncollinear configurations are stable. [2ll. [2^) The equa- 
tions Ji > 0, Ji > 0, Ji — 2Ji > determine the region in the parameter space where a 
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FIG. 1: Basic spin-configurations in isotropic antiferromagnetically coupled two-layers (A^=2). At 
zero field (a) the magnetization vectors rrij are antiparallel in adjacent layers but the structure 
has no fixed orientation in the plane (infinite degeneracy). An in-plane magnetic field lifts the 
degeneracy and stabilizes states with 1 perpendicular to the field and with m along the field - this 
is a spin- flop (SF) phase (b). In a sufficiently strong magnetic field the SF phases transform into 
the saturated (flip) phase (c) 

coUinear antiferromagnetic (AF) phase is the zero-field ground-state consisting of blocks of 
adjacent layers with antiparallel orientation of the magnetization (Fig. 1(a)). We restrict 
our analysis to this practically important case. The AF phase is infinitely degenerate with 
respect to orientations of the staggered vector 1 in planar directions. 

The magnetic states for the isotropic two-layer (A^ = 2) in the applied field can be 
obtained by minimization of the energy ^ with K = 0. The solutions {6 = ip, H = 
2(J — 2 J) cos + 8Jcos^0) describe the states with the staggered vector perpendicular 
to the applied field (so called spin-flop (SF) phase). In an increasing magnetic field the 
magnetization vectors rotate into the field direction (Fig. 1 (b)), and finally a continuous 
transition into the ferromagnetic {flip) phase with ^ = 0. This {spin-flip transition) occurs 
at the exchange-field He = 2{J + 2J) (Fig. 1 (c)), which is lower than the corresponding 
spin-flip field = 4( J + 2 J) for a bulk system. 

In superlattices with A^ > 4 the magnetic field applied in the plane also lifts the degen- 
eracy and stabilizes the SF phase (the state with 1 ± H and 1||H) (Figs. 2 to 4). However, 
the magnetic configurations in these SF states and their evolution in the applied field are 
markedly different from those in bulk antiferromagnets or in the two-layer systems. It turns 
out that magnetic structures of the multilayers with A^ divisible by four (A^ = 4/ called 
here even-even) differ from even-odd systems with A^ = 4Z + 2 (/ = 1,2...). In low fields ( 
Ji — 2 Jj ^ H) the solutions for the SF phase are given by a set of linear equations 

{J2j-i-2J2j-i){n-e2,-, + e2j) = H, e2,-e2,+i = o, j = i,2...z, (4) 

where I = N/4 for even-even and I = {N + 2)/4 even-odd systems. It is clear that for the 
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isotropic model the direction of the magnetic field in the film plane plays no role. Following 
the definition of the angles 6i we assume here that the magnetic field is applied along the 
X-axis. The solutions (jlj) describe small deviations of the magnetization vectors, from the 
directions perpendicular to the easy axis (see the magnetization profiles in Figs. 2, 4 and the 
configuration in Fig. 3 (b)). In all internal layers {i = 2....N — 1) neighbouring pairs retain 
antiparallel orientations (e.g., (m2,m3) and (1114,1x15) in Fig. 3b). For even-even systems 
the magnetization vectors of the central layers (mjv/2) '^N/2+1) preserve perpendicular ori- 
entation even in finite (but weak) fields (see layers 6 and 7 in Fig. 3). For even-odd systems 
the magnetization vectors of all layers have finite deviations from perpendicular orientation. 
Towards top and bottom layer i = 1 or in the stack, these deviations increase. Due 
to antiparallel orientation the pairs have zero net magnetization. The total magnetization 
arises only due to the top and bottom layers i = 1 and A^. Note that in these configurations 
a number of layers rotates against the applied field (in Fig. 2-4 the layers with z = 2, 4 
and their symmetric counterparts with i = 11, 9). This occurs because, in weak fields, the 
exchange interactions favouring antiparallel magnetizations in adjacent layers play the dom- 
inant role. An increasing magnetic field counteracts and slows down this reverse rotation, 
when the non-linear evolution for the magnetization structure is reached, and finally the 
sense of rotation is changed at characteristic fields Hr'^ where dOi/dH = 0. With increasing 
A^ the number of layers increases which display this reverse rotation, the deviations from 
orientations perpendicular to the applied low field become larger near the surface layers, and 
the fields H^'^ are reached at lower fields (see inset in Fig. 2 for A^ = 64). A further set of 

characteristic fields H^f^ defines the points where the projection of nij on the field directions 

(2) (2) 

changes the sign (the field Hr and H]_ are indicated in Fig. 2.) In increasing field these 
characteristic fields initially are reached for the central layers and at higher fields for those 
closer to the boundaries. For fields H > Hf^ = H^^~^^ the magnetization directions of all 
layers have positive projections onto the field direction. In the model with equal exchange 
constant there is another special field H* (independent of A^), where all inner layers have the 
same projection onto the field direction {6i = (— 1)*+^^^Q) * = 2,3...A^ — 1). The parameters 
of this "knot" point are determined from a set of equations 




K = 2J/J. 
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FIG. 2: Evolution of spin-states in multilayers with a field H applied in the plane of the layers 
and for zero anisotropy K = Q and biquadratic exchange J = 0, hence = 4 J. Main figure: 
rotation angles of the layers in the plane 6i i = 1, . . . N for super lattice with = 12. At point 
labeled "R" the sense of rotation changes for the magnetization direction 9i (here, i = 2); at point 
"_L" the magnetization of this layer is perpendicular to the applied field again. Point 'V is the 
"knot-point" (Eq. @). Inset: = 64 only 9i for even i = 2 and i = 4, 8, 12 ... 64 shown. 

The functions H*{k), 01{k,), 6*{k) are plotted in Fig. |S| In particular for superlattices 
with zero biquadratic exchange J = one has H* = V6He/4: and 9* = acos (v/378). Near 
saturation, the SF phase has only positive projections of the magnetization on the direction 
of the magnetic field which decreases towards the center similarly to spin configurations 
obtained numerically in Ref. 13 

The spatial inhomogeneity of the SF phases in the multilayers and their remarkable evo- 
lution in the applied field are due to the particular finite-size effects in this type of magnetic 
nanostructures. In such layered nanostructures all internal layers interact with two adjacent 
layers but the top and bottom layer have only one neighbouring layer. Thus, their ex- 
change coupling is weakened compared to internal layers within the multilayer stack. This 
strong disbalance of the exchange forces at the boundaries affects the magnetic ordering 
within the entire multilayer and causes the reorientational processes in the SF phases. This 
mechanism is also responsible for the occurrence of inhomogeneous states in other non- 
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FIG. 3: Spin-states in multilayers with = 12 for zero anisotropy K = and field H applied 
in the plane: (a) H = (b) H < ifjj^^ <C He (c) H = (d) H = Hf^ (e) H = H* , and 

(f) H* < H < He'^\ i.e. below the transition to the flip-phase (cf. Fig. |2j. All states have 
mirror-symmetry about the center of the multilayer stack (marked by the dashed line). 

collinear structures in the multilayers (see IV. C and V.A). Note, this important role of the 
cut exchange-couplings is specific to antiferromagnetic confined systems with noncompen- 
sated moments at surfaces. In contrast, the cutting of exchange couplings does not affect 
ferromagnetically coupled nanostructures where relative orientation of the magnetization in 
the layers essentially depend only on magnetic anisotropy, applied fields, and the demagne- 
tizing stray-fields. 

B. Effects of the biquadratic exchange 

Evolution of spin-flop profiles for the multilayer with finite biquadratic exchange is plot- 
ted in Fig. 4. For the systems with collinear ground-states the biquadratic exchange does 
not induce reorientional transitions but it has rather strong influence on the distribution of 
the magnetization in the multilayers (see insets in Fig. 4) and the value of the character- 
istic fields (Fig. 5). The effects of a biquadratic exchange J > 0, if present in multilayers 
with collinear antiferromagnetic ground-state, can be easily understood. For magnetic con- 



12 




FIG. 4: Evolution of magnetic states in antiferromagnetic multilayer = 12 with zero anisotropy 
K = and biquadratic exchange n = 2J/J > 0, i.e. with enhanced Hg = 4(J + 2J) (compare with 
Fig. [21 for «; = 0). H is applied in the layer plane. Insets (a): magnetization curves in the range of 
values of k with antiferromagnetic ground-state in zero field, (b): evolution of states (only shown 
for layers with even numbers) in small fields for k = 0.98. 

figurations close to the antiparallel orientation of neighbouring layers, it softens the linear 
exchange forces; for nearly parallel orientation the system becomes stiffer instead. The evo- 
lution of the magnetic states in the external field is distorted accordingly by the presence 
of the biquadratic exchange J > (Fig. In low fields close to the antiferromagnetic 
ground-state, the system reacts more strongly on the external field. Thus, the fields Hr\ 
where the rotation of the i'the layer reverses, and the field H^l\ where perpendicular ori- 
entation with respect to field is regained, respectively, are reached at relatively lower fields 
(labels "R" and "_L" in Fig. 0}. Conversely, the ferromagnetic state is reached only close to 
the enhanced exchange fields He reflecting the stiffening of the effective exchange couplings 
in nearly ferromagnetic configurations. Thus, the magnetization curves acquire a nonlinear 
character with increasing J > (inset (a) in Fig. E}. Experimental observation of the fields 
Hr^ and the special field H* (Fig. could be used to determine the relative strength of 
J > 0. For large values of J the reversal of the rotation of the magnetizations of individual 
layers takes place at very low fields (see, inset (b) in Fig. H}. 
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FIG. 5: Dependence of "knot" point H* (Eq. © on strength of biquadratic exchange k. Inset 
shows the two angles 9*, 9\ which characterise this special configuration independent on N . 
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FIG. 6: Spin- flip fields for antiferromagnetic multilayers (zero anisotropy = 0) in dependence 
on number of layers A^. 

C. Spin-flip transition 

Near the spin-flip transition from the SF-phase to the ferromagnetic states, the deviations 
of irij from the field directions are small {9i <^ 1), and the energy of the system can be 
expressed by the quadratic form W = X]f^=i ^ij^i^jy where the matrix Aij has a tridiagonal 
band structure with nonzero elements only in the main diagonal An = H — Ji and with the 
side diagonal elements = Ai^i^i = Ji + 2Ji = Ji. The spin-flip field ifi^'* is determined 

as the largest root of the equation det{Aij) =0. In particular for the model with equal 
parameters these solutions are plotted in Fig. 6. The spin-flip fleld gradually increases from 
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hP = 2(J + 2J) for the two-layer to the "bulk" value He = 4(J+2J) in the limit of infinite 
(Fig. 6). This dependence refiects the increase of the average number of exchange bonds 
from the value 1 for the two-layer and approaching 2 as iV tends to infinity. 



IV. EFFECTS OF TETRAGONAL ANISOTROPY IN TWO-LAYER SYSTEMS 



Effects of the four-fold magnetic anisotropy on the states in the antiferromagnetic super- 
lattices are revealed in detail by applying fields in arbitrary directions. We approach the 
general case by a detailed investigation on the case of the two-layer. To set the stage, we 
discuss the highly symmetric phase-diagrams with fields in direction of easy and hard-axes 



in the plane. This completes earlier work by Dieny et al. [2j| on the N=2 systems. Then, we 
present phase diagrams for arbitrary field directions in the layer plane. As we have seen in 
the last section, the primary effect of biquadratic exchange is a distortion of phase diagrams 
in various regions. Thus, here we consider only models with J = to avoid such quantitative 
effects which will not affect the general topological features of the phase diagrams. 

Optimization of the function $2(^5 0) Eq- © yields solutions for the magnetic states and 
their stability limits. In zero field the tetragonal anisotropy lifts the rotational degeneracy 
of the AF phase and stabilizes two different states with perpendicular orientations of 1 (Fig. 
7 (a)). 



A. Magnetic field along hard directions 

In a field applied along one of the hard axes, the AF configurations of the ground-state 
(Fig. 7(a), (b)) are distorted into low symmetry configurations (so called canted phases - 
Fig. 7(c), (d)). Both these magnetic configurations preserve mirror symmetry with respect 
to the field direction and have the same energy. The total magnetizations m of these canted 
states deviate from the field directions to different sides at equal angles (Fig. 7(c'),(d')). 
An oblique magnetic field deviating from the hard direction to one or the other side violates 
the phase balance between these two canted states. It favours the canted phase with the 
larger projection of m onto the field direction. This is the typical situation of a first-order 
(or discontinuous) transition. 

In an increasing field the vectors m in both canted phases rotate into the field direction 
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FIG. 7: Magnetic configurations in planar antiferromagnetically coupled two-layers with four- 
fold anisotropy. Easy directions are indicated by thin continuous lines, hard directions by dashed 
lines. In zero field two antiferromagnetic configurations (I) (a) and (II) (b) with the staggered 
vectors along the easy directions correspond to the ground-state. A magnetic field along the hard 
directions distorts them into the canted phases (c), (d) which transform into the spin-flop state 
(e). Corresponding configurations of net magnetization and staggered vectors are shown for the 
two canted phases (c') (d'), and for the spin-fiop phase (e'). 

and reach it at a certain critical field Hs/i- At this point both phases transform into the 
common configuration corresponding to the SF phase (Fig. 7(d)), i.e. a second-order (or 
continuous) phase transition from low symmetry canted phases into the high-symmetry 
SF-phase occurs in this field. Standard analysis yields the following expression for the 
parameters of the critical point 

H,f^ = {J + Kiy)^/2{l + iy), u = {1 + k - Vl + Uk + 25P)/(6A;), k=\K\/J, (6) 

where u = cos (psf determines the equilibrium SF configuration at the critical field. 

The two canted phases are competing phases related by a first-order transition. At the 
field Hsfi the discontinuity between these phases disappears, i.e. this is a critical end-point 
of the first-order transition. After the transition into the SF phase, the system further 
evolves by rotation of the magnetizations mj into the direction of the field. This process is 
terminated at another critical field Hfi = 2J(1 + k) (Fig. 8). 

The transition between the canted and SF phases is continuous only below a certain 
anisotropy strength, \K\ < K^, and becomes discontinuous at higher anisotropy. The 
parameters of the corresponding tricritical point Ka = 0.081406 J, Ha = 0.425780 J have 
been calculated numerically from the equation = (1 + z/)/ (6z/^ — 2z/ — 4) together with 
Eqs. (P). Numerically calculated transition fields Htri and the upper stability fields of the 
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FIG. 8: Magnetic phase diagram of the two-layer system in magnetic fields along hard directions. 
The tricritical point {Kti = 0.081406 J, Hti = 0.425780 J) separates the continuous (dotted line) 
and discontinuous (thick line) transitions between the canted and SF phases. A further dotted line 
gives the critical field of the second-order transition between the SF and flip phases Hfi Thin 
lines give the lower stability limit H^fi of the SF phases and the upper limit Hd of the canted 
phases. The grey area is the region of the phase equilibrium between the two canted phases (see 
Fig. [3 (c),(c') and (d),(d')). Magnetization curves for low and high anisotropy, as indicated in the 
main panel, are plotted in (a) and (b). 

canted phases Hd are plotted in Fig. 8. In this picture the grey area describes the region 
of the first-order transition, respectively the region of phase-coexistence between the canted 
phases. 

B. Magnetic field along easy directions 

A magnetic field applied in one of the easy directions violates the energy balance between 
the two AF states. The state AF(I) with staggered vector parallel to the field does not change 
its configuration and exists as a metastable collinear state up to the critical field Haf. The 
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phase AF(II) with the staggered vector perpendicular to the magnetic field transforms into 
a SF phase which corresponds to the global energy minimum of the system. In an increasing 
magnetic field this phase evolves similarly to the SF phase along the hard directions, and 
continuously transforms into the fiip states at the critical fields Hf2 for low anisotropy. 

The transition between the SF and fiip phases changes the order in another tricritical 
point {Kt2,Ht2)- At a critical field Hsf2 © the SF phase becomes unstable with respect 
to fiipping distortions, i.e. modes that redirect the magnetization vectors into the field 
direction but preserve the symmetry of the SF phase. The parameters of the tricritical 
point (|-ft'|t2 = (1/5) J,Ht2 = (8/5) J), and the stability limits 



H f 



2Jv/Ml+^, ^^/2 = ¥V^^^^' Hf2 = 2J{l-k) (7) 



have been calculated in |22i]. After some algebra the first-order transition field Htr2 also can 
be derived in analytical form 



Htr2 - - 




(8) 



At the transitional field we have cos = (^1 + 3/A;- l)/3. The (K, H)-phase dia gram for 
this cases is plotted in Fig. 9. 

With increasing strength of the anisotropy, the "landscape" described by the energy 
function <1>2(6',0) acquires additional folds. In particular, metastable canted phases arise in 
a certain region of the magnetic fields along the easy directions. (This region is indicated by 
grey colour in Fig. El). For \K\/ J = k > k* with k* = 0.28+0.08 Vg = 0.47596 the SF phase 
undergoes an instability with respect to a transition into the canted phases. This instability 
occurs at lability fields which are given by the two branches of the following parametric 
equation 



ifijf = (J - \K\u,^2) ^2(1 + ^1,2), ^1,2 = {k-lT Vl-Uk + 25k^)/{6k), (9) 

where i>i 2 is given by the the configuration in the SF-state, z/i 2 = cos20sj2- Both branches 
of the lability field © start at the point H^^ff{k*) = 1.3895 J (point □ in Fig. ij; H^f^ 
meets the lability line Hsf2 of Eq. © at the "beak" {k = 1/2, H = V2 J) (point A in Fig. 

lO}. The transition between the metastable SF phase and the metastable canted phases along 

(I 2) 

Hlj2 is continuous between the points A and a further tricritical point Q- For stronger 
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FIG. 9: Magnetic phase diagram of the two-layer system in magnetic fields along easy directions. 
The tricritical point [Kt2 = (1/5) J, Ht2 = (8/5) J) separates the continuous Hf2 (dashed 
line) and discontinuous Hf2 ©(thick line) transitions between SF and flip phases. The tricritical 
point [Kt2 = 0.2J, Ht2 = 1.6J) separates the continuous Hf2 ^ (dashed line) and discontinuous 
Hf2 ©(thick line) transitions between SF and flip phases. Thin lines indicate the lability lines 
of the individual phases Q. The grey area designates the stability region of the canted phases. 
For the anisotropy \K\/ J = k* = 0.47596 (point □), the critical field of the SF phase switches 
from an instability against symmetric "flipping" at the line to an instability against canting 
at given by the two branches Hgj^ from Eq. Q. (Inset in main figure shows the details of 
this process; see text following Eq. ® for further explanation.) Small panels (a) and (b) show 
magnetization curves for low and high anisotropics as indicated in the main figure. Full thick lines 
are for the evolution of the equilibrium ground- states, in (a) from SF phase to a continuous spin- flip 
transition into saturation. Dotted lines are the evolution starting from the metastable AF-phase. 
In (b) the evolution of the saturated state in decreasing fields is given by a thin line. The inner 
small hysteresis loop shows the transitions from and into the metastable canted phase. 
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anisotropy the transition between these metastable phases is discontinuous at a first-order 
transition field Htr^- This transition fine along with the stability limits of the canted phases 
have been numerically determined as given in Fig. El Note that such processes between 
metastable phases may be realized only in systems with high coercivity where the transition 
into the thermodynamically stable fiip phase is suppressed at the field HtT2 (0)- 



C. Evolution of {H^, Hy)--phase diagrams 

Now we consider the magnetic states in applied fields deviating from the symmetric 
directions and construct (if^, ify)-phase diagrams for different values of the parameter k = 
\K\/J. For the (if^;, ifj,)-diagrams in Figs. ITUl to 14 it is assumed that the easy directions 
coincide with the x- and y-axis. This corresponds > in the energy Eq. Q. 

Diagrams for K < can be obtained by rotation of those for > through an angle 
7r/4. In the limit of weak anisotropy, \K\ <^ J, independent minimization with respect to 
the angle (p yields 2 J cos = H cos{6 — ip) and the energy (jHl) is simplified to the following 
form 



F{e) = cos2(^ -^)-— COS 4^ 



2 



l + ^cos^(^-V^) 



(10) 



The set of equations for the stability limits of solutions, dF{6)/d6 = and (PF{6)/d6'^ = 0, 
yield two closed lines of critical fields in the {Hx, ify)-phase diagram (Fig. EH)- One of them 
H = 2{J — K cos Alp) describes the second-order transition into the fiip phase with = 0. 
The other closed curve can be written in a parametric form 

cos(2^ - 2^) = -AKJ cos AO + K^fli (9) , 

sm{2e - 211;) = -2KJ sinAO + K^VL2{e) , (11) 



where fii(0) = AK"^ {7 cos^ AO - 2 + 2sgn(K) VScos^i^ + l, fi2(^) = 2K'^{bcosAe + 



2sgn(i^') V3 cos^ 46' + 1. It describes an astroid with eight-cusps that confines the region 
of the canted metastable states (Fig. imT a)). The cusps along the easy directions coincide 
with the stability field of the AF phase Haf Eq. 0, and those along the hard directions 
with Hsfi Eq. Q. Thick black lines within the astroid and along the hard directions in- 
dicate the first-order transitions between the canted phases from Fig. [7| This topology of 
the (ifj,, i/y) -phase diagram Fig. irUT a) is preserved up to the first tricritical point Kn- For 
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FIG. 10: (a) Schematic {Hx, Hy)-phase diagram of the two-layer for < K < Kti- The astroid 
(continuous hne) gives the stabihty Umit of the canted metastable phases. Thick hnes in diagonal 
directions (hard axis) are first order transitions between canted phases. The line for the (continu- 
ous) transition into the saturated spin-flip state is given by the dotted line. Note that the real size 
of the astroid in the {H^, Hy)-dia,gram is much smaller than shown here, (b) Evolution of spin-flip 
field, respectively lower stability limits of the saturated (ferromagnetic) phase (FM). Phase FM is 
stable outside the closed curves for the various values of anisotropy < K < 1. For K/J > 1.0 
the FM-phases are (meta)stable even in reverse fields and the existence regions overlap for states 
saturated in different easy directions. 

I-^I > a first-order transition between SF phase and two canted phases arises along 
the hard directions (point a in Fig. ITTT a)). For stronger anisotropics, \K\ > Kfi, discon- 
tinuous phase transitions exist also for finite deviations of the field direction from the hard 
axes. The corresponding line of these first-order transitions have critical end-points, where 
the difference between the competing phases disappears (analogously to the end-point of a 
coexistence-line in gas-liquid phase diagrams). For \K\ > Kt2 and for fields along the easy 
directions, the spin-flip into the saturated (induced) ferromagnetic state occurs discontin- 
uously. Thus, another line of first-order transitions develops for \K\ > Kt2 also in oblique 
fields. Correspondingly, the transition lines for spin-flips in (if^., ify)-phase diagrams consist 
of continuous and discontinuous sections joined by the tricritical points 5, 5' . The calculated 
(i^a;, i7y)-phase diagrams for K = 0.3 J in Fig. [TTl (a) and for K = 0.8 J in Fig. [TT] (b) rep- 
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FIG. 11: First quadrant of the {H^, Hy)-phase diagram for antiferromagnetic two-layers with (a) 
K = 0.3 J (b) K = 0.8 J. Easy axes of magnetic anisotropy are along the axes of the diagrams. 
Inset in (a) magnifies the region of the "swallow tail" for the discontinuous transitions in magnetic 
fields oriented close to the hard axes. Thick lines are lines of first-order transitions. The dotted 
lines 5-5' are continuous spin-flips from spin-flop phase to the saturated (spin-flip) state. Thinner 
lines are the limits of stability, deflning the coexistence regions for the corresponding phases. 

resent the main topological features of these phase diagrams. Point 'a' indicates the field 
of phase equilibrium between the SF phase and two canted phases and corresponds to the 
critical field Hfri from Fig. 8; Points (3, (3' are the end-points of the first-order transitions 
between one of the canted phases and the SF phase. Points 7 7' correspond to the transi- 
tion field Htr2 from Fig. 9; points 5, 5' are the tricritical points, where the discontinuous 
transition between the flip phase and the distorted SF phase ends. The lability lines in 
the vicinity of the transition field Htri have the shape of a "swallow tail". Similar phase 
diagrams arise in uniaxial ferromagnets with strong fourth-order anisotropy. 0, Q As the 
anisotropy constant increases in the region l-ft'l > Kt2 the lines of the first-order transitions 
extend and near the special value \K\/ J=k* = 0.476 the points f3, (3' reach the easy direc- 
tions. For larger \K\ canted phases exist as metastable state also for fields in easy directions 
(compare Fig. 9). However, for increasing fields close to this direction, the discontinuous 
transition from the spin-flop state into the spin-flip state occurs before the transition into 
such a canted phase can take place, i.e. the two first-order lines a - (3 and 5 - 7 cross each 
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other close to the easy axis directions in Fig. 11(b). Brohl et ah [0\ reported evidence of an 



d in an easy direction, 
2211 However, the state 



intermediate canted state in a Co/Cu/Co/Cu(001) system with fie^ 
following the theoretical prediction of such states by Dieny et al. 
was found at lower fields than expected. This may be due to a misorientation of the external 
field and/or a mosaic of the epitaxial layer system because then the canted state becomes 
stable already at lower fields, as can be seen from line a - /? in phase diagram Fig. 11(b). 
Generally, it is not sufficient to investigate only magnetization behaviour in hard and easy 
axes directions for a thorough understanding of the magnetization phases in these multilayer 
systems. 

To understand the transformation of {H^, ifj^)-diagrams with increasing tetragonal 
anisotropy the limiting case is useful where rrij are strictly oriented along the easy axes. 
In this case with infinite anisotropy, i.e. \K\ = oo, our model can be considered as a chain 

in a transverse 

external field. Recently, such four-state-clock models were employed to analyse the spin 
configurations of bulk tetragonal metamagnets with large four-fold anisotropies, such as 



rare-earth nickel-borocarbides 







4l| or rare-earth silver-antimonides. 42] Rich experimen- 



tal (if^., if y) -phase diagrams have been obtained and analysed in terms of four-state-clock 
models, e.g. for HoNi2B2C [41j or DvAgSb9.|3| The values of the angles 6i = 0, tt/2, tt, or 
3 7r/2 in these four- fold states for i^' = oo are symbolically given by f, |, and e.g. the 
AF and SF phases are (t |) and {—>■ ^) for N = 2. In addition to the coUinear AF and flip 
states there are phases with perpendicular orientation of the magnetization, or "90 degree- 
folded" phases (Fig. I12|l. The states created by all combinations of these four angles 6i for 
the magnetization in the multilayer stack exist as metastable states in arbitrary magnetic 
fields because the energy wells corresponding to these solutions are separated by infinitely 
high potential barriers. The regions of the absolute stability of these phases are separated 
by first-order transitions lines shown in Fig. El In the first quadrant the transition between 
the AF states and the canted phase (t —>■) occurs at the line + Hy = J; the transitions 
lines from the canted phase (t — ;>) into the flip states (t t) and (— — *>) are Hy — = J 
and Hy — H^ = —J, correspondingly (Fig. 12). In the points (1,0) and (0,1) four phases 
coexist. Thus, under increasing strength of the anisotropy {K > 0) the {H^, Hy) -phase 
diagram evolves from that plotted in Fig. [TUT a) for K <^ J to that in Fig. for infinitely 
large values of K. 
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FIG. 12: {Hj., Hy)-phase diagram of the N = 2-layer system in the hmit of an infinite positive 
four-fold anisotropy. 

V. MAGNETIZATION CONFIGURATIONS AND PROCESSES IN MULTILAY- 
ERS 

The equihbrium magnetic configurations in two-layer systems arise as a results of the com- 
petition between the interlayer exchange coupling and tetragonal anisotropy. For multilayers 
the disbalance of the exchange forces at the boundaries (Sec. II) additionally influences the 
magnetic states. We first describe the structure of the one-dimensional solutions for lat- 
erally homogeneous states in finite antiferromagnetic superlattices. These magnetic states 
are determined by the interplay between cut exchange at the surfaces and the four-fold 
anisotropy. Then we discuss some consequences for multidomain states and magnetization 
processes, and we discuss the physical nature of other effects which may play a role for the 
magnetic behaviour of real experimental systems. 

A. Exchange cut versus tetragonal anisotropy 

First, we investigate the {H^, iJy) -diagrams of magnetic states of multilayers with N > A 
in the limit of infinite four-fold anisotropy. The zero-field ground-states are the AF phases 
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FIG. 13: {Hx,Hy)-phase diagram of the multilayer in the limit of infinite positive four-fold 
anisotropy. The phase-region of the L'^^^X'^^^-phases coincides with that of the collinear ferri- 
magnetic F^-^^, not shown here. (For details see text.) 

(T i ■ ■ ■ T i) and <— ■ ■ ■ ^ ) We may restrict the field to be oriented in directions 
ijj in the interval [0,7r/4]. The configurations for other values of ijj follow from symmetry. 
Then, we may distinguish the four ground-state (zero- field) domains (AF1)(")= (Ti)^"^ 

(AF2)(")= (it)^"\ (SF1)(")= (^^)("\ and (SF2)(")= (< ^)("), where n signifies the 

number of repetitions of the pair in a domain. In external fields H > configurations 
with net magnetization will be stabilized. These configurations must have pairs of adja- 
cent moments flipped by 90 degree L=(t — >) or X=(<— f), and by 180 degree F=(t t)- The 
90-degree-folded configuration "X" is less favourable than "L" for field orientations in the 
chosen range, except for fields in easy- axis direction ip = 0. Thus, the magnetized con- 
figurations with lowest magnetization and smallest expense of exchange energy arc those 
with only one pair of type "L". In a short-hand, we write L^^^ for these configurations of 
type (AF1)(") /L /(SF2)(^/2-'*-i) (with n = 0, 1, . . . .) Next, we may form configurations 
with higher magnetization and smallest expense of exchange energy using one "L" and one 
"X"-pair: L«X«= (AF1)(") /L /(SF2) W^-"-"^-^) /X /(AF1)("*), (n,m = 0, 1, . . .). These 
configurations have the same energy and, therefore, are degenerate with configurations con- 
taining only one "F"-pair: F(i) = (AF1)("VF /(AF2) W^-"-^). Note that these states are 
highly degenerate because the "L"-, "X"-pair, or "F"-pair may be placed at arbitrary posi- 
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tions in the stack of layers. In higher fields configurations with various combinations of 
"L"-, "X"-, and "F" -pairs may be stabilized by an external field. However, the structures 
with lowest energy, i.e. the absolutely stable states, are rather simple because of the fol- 
lowing considerations. If, in external fields, more than one "L"-pair can be formed starting 
from the L*^^) structure, the formation of the maximum number of L-pairs gives the most 
favourable configuration. In particular, the structure L'^^/^) = (| — is the lowest en- 
ergy structure for fields pointing in hard axis-directions ip = n/A in the limit of infinite H. 
The state with higher magnetization than L*^^/^^ for fields pointing closer to the easy axis 
direction, < ip < it /A, and with the smallest expense of exchange energy is in our nota- 
tion L'^^^'^^^'^F. The saturated state FM=F*^^/^) is the most favourable state which is reached 
whenever energy can be gained in external fields by fiipping more moments into position | in 
the states L^^/^~^^F, F'^^^ or L^^^X^^^^. Thus, for general N we have only the following phases 
for infinite positive four-fold anisotropy: AF or SF as zero-field ground-states; for fields with 
orientation close to the hard direction ijj = n/A L^^-* and L*-^/^-*; otherwise in intermediate 
fields two degenerate phases F^^^^ and L'^^^X^^^; at high fields an asymmetric L*^^/^~^)F-phase 
and the fully saturated ferromagnetic phase FM. We did numerical checks to ascertain that 
no further energetically stable phases do exist in external fields of arbitrary strength and 
direction, indeed. Thus, we searched for the states of lowest energy by sampling all possible 
configurations for models with N=4, . . . , 12 corroborating our arguments. Based on this set 
of magnetic configurations the resulting {Hx, Hy)-ph.eise diagram for general > 4 can be 
calculated analytically (Fig. fT!?|). As in the simpler case of Fig. El for the two-layer system 
all these states are separated by infinitely high potential barriers and remain metastable 
for arbitrary fields. The first-order transitions between different phases occur along straight 
lines as shown in Fig. 

For finite strength of the four-fold anisotropy and under the infiuence of the exchange 
interactions, the basic structures are derived from the phases in Fig. El Under the infiuence 
of the field, they are elastically distorted into spatially inhomogeneous configurations. We 
have numerically investigated models for such cases with A^ = 4 ... 20 and various values of 
anisotropy K. Fig. El displays the general features for the example A^ = 8. The numerically 
calculated magnetization curves corresponding to the lowest energy states for sufficiently 
high anisotropies show the sequence of phases present in the infinite-anisotropy phase di- 
agram Fig. 1131 However, the degeneracy of these phases is lifted because the distortions 
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possible at finite anisotropy yield different gains of energy for the different configurations of 
the phases L^^^^ F'^^\ and L^^^X'^^^ Generally, depending on field-orientation, large domains 
with nearly spin-flop-like configurations 1 ± H are favoured because these configurations can 
be more easily distorted by the fields yielding a corresponding gain of energy. In particular, 
phases derived from the ferrimagnetic collinear F'^^-'-phase are disfavoured compared to the 
L(i)X(^)-type phases. The nearly collinear ferrimagnetic configurations seem to exist only 
as metastable states. For fields closer to the hard-axis direction the phases with either one 
L-pair or N/2 L-pairs occur. Under the influence of an applied field further discontinuous 
transitions occur at intermediate anisotropics, as seen e.g. in Fig. El (a) for K/ J = 0.375. 
These transitions are jumps from one energy basin to another which are formed by distorting 
degenerate configurations mainly of the type L^^^^ and L^^^X^^^ 

At low anisotropics the phase diagram attains the behaviour discussed for zero 
anisotropics in section III. Here the magnetic states are only influenced by the cuts of 
exchange bonds at the surfaces and spatially inhomogeneous spin-flop states slightly dis- 
torted by four-fold anisotropy are realized (see Figs. 2 to 4). (if^;, iJj,)-diagrams for this 
isotropic case consist of the region of the inhomogeneous SF phase (Fig. 2, 3) separated 
from the saturated flip state by the critical line He^^ (Fig- 6). These global features of 
the phase-diagrams are similar to the case of antiferromagnetic multilayers with uniaxial 
anisotropy. ^ 

Up to now we have discussed the effects of the competition between the bilinear exchange 
and four-fold anisotropy. It is clear that finite biquadratic coupling may substantially change 
the critical fields and stability regions of the various phases as it affects the elastic stiffness 
of the system. We note further, that the simple structure of the phase-diagram ruled by 
the phases present in the infinite anisotropy limit of Fig. is valid only for the case of 
equal exchange constants Jj = J in the multilayer stack. For arbitrary sets of values for Jj 
in energy the phase diagrams may become considerably more complex and may contain 
further phases with different combinations of flipped 90 and 180 degree spin-pairs. Even 
in such cases, the outline of the magnetic phase-diagrams, described here for finite equal- 
constant superlattices, should generally hold: (i) The high anisotropy limit is comparably 
simple with few phases determined by the competing lowest-energy basins of the anisotropy. 
Non-equal constants may lift some degeneracies that are present in superlattices. (ii) In 
intermediate anisotropy range many elastically distorted phases appear, which are derived 
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FIG. 14: Examples of magnetization in an antiferromagnetic multilayer with N = 8. Curves 
corresponding to the evolution of lowest energy states and for various values of anisotropy are 
shown with field oriented in directions of an easy axis (a), in an oblique field (b) and in hard-axis 
direction (c). ip is angle between magnetic field and easy axis direction. In each case, various 
magnetic configurations are shown for the discontinuous evolution of the magnetization curves 
with highest anisotropy. 

from stable and metastable high-anisotropy phases, (iii) For vanishing anisotropy the phase- 
diagrams become simple again as only inhomogeneous spin- flop-like phases and the saturated 
ferromagnetic phase remain in external fields. 
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B. Magnetization processes in real systems 



So far, we have analyzed single domain magnetic configurations. Magnetization pro- 
cesses in real systems, however, are usually accompanied by complex reconstruction of the 
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multidomain patterns as those recently observed in Fe/Cr multilayers. 

There are two main physical mechanism for multidomain states in the systems under 
discussion. Four-fold degeneracy of the ground-state leads to creation of antiferromagnetic 
multidomain structures with 90 degree domain walls. Unlike the case of magnetic materials 
with nonzero total magnetization, where multidomain states are caused by demagnetization 
effects, in antiferromagnets such domains are metastatic and arise during the formation of 
the ordering state, i.e. they have a kinetic origin. Hence, multidomain patterns observed 
in antiferromagnetic coupled two-layers and multilayers have irregular morphologies and 
depend on thermal and magnetic- field histories, j^] 

However, another type of multidomain structures arises in the vicinity of field-induced 
discontinuous transitions, Such thermodynamically stable transitional domain structures 
are formed by domains from states corresponding to the coexisting phases at first-order 
transitions. These domains are analogous to the domains of a demagnetized ferromagnet. 
In principle, the equilibrium parameters of such multidomain structures and their boundaries 
can be calculated by standard methods. Q,Q In Q such calculations have been carried 
out for bulk easy-plane tetragonal antiferromagnets. Magnetoelastic interactions lead to 
modification of the inhomogeneous states and decrease the regions of the multidomain states 



up to their complete suppression 



ion.Q, 
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44l | This and other coercivity mechanism partly 



block the development of the equilibrium states. As a result, in real systems the evolution of 

ather strong hystereses. j9(] For experiments. 



multidomain states 



there are two important consequences related to the starting states and history dependence 
of magnetization processes. The "texture" of a real antiferromagnetic state in zero-field 
depends on the detailed kinetics imposed by e.g. cooling rates or deposition conditions. 
On the other hand, field cycling by inner loops for an antiferromagnetic multilayer may 
yield various metastable configurations and multidomain structures owing to the very wide 
coexistence region of the domains whenever sizable magnetic anisotropies are present in the 
multilayer stack. 
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VI. CONCLUSIONS 



Within a phenomenological approach we develop the theory of reorientation transition in 
antiferromagnetically coupled superlattices with in-plane magnetization. Detailed investi- 
gations of the surface effects in the isotropic multilayers (Sec. Ill) and four-fold anisotropy 
effects in two-layer systems (Sec. IV) reveal the most important features of the system: (i) 
Complex evolution of the inhomogeneous states (Figs. 2-4) is imposed by the strong dis- 
balance of the exchange coupling, i.e. by the cut of the exchange bonds, (ii) Remarkable 
field-induced reorient ational processes occur due to enhanced four-fold anisotropy (Figs. 8 
to El). 

The model used here, corresponding solutions and phase-diagrams include as special 
cases earlier theoretical studies on surface j3] and four-fold anisotropy 



22, 29^ effects. The 



results are in essential accordance with existing experimental observations on inhomogeneous 
distribution in the spin-flop phase near the saturation field and some effects of four- fold 
anisotropy in Co/Cu(001) wedged two-layers ^ and Fe/Cr superlattices. [3] Our approach 
and results enable a qualitative analysis of the magnetization processes in the multilayered 
systems (Sec. V). In spite of the rather complex phase- diagrams of these systems, the 
analysis can be extended towards a quantitative description of real systems belonging to the 
class of artificial layered antiferromagnets described by Eq. (^. 

So far most experimental work is carried out only for special conditions, often data are 
collected only with fields along easy-axes. These results only cover small regions of the 
{Hx,Hy)-ph.a.se diagrams (Figs. ITUl to IT^ and do not capture the rich varieties of magne- 
tizations processes in such systems. It is remarkable that many interesting effects , such 
as 90-degree folded phases, transitions into asymmetric canted phases etc., are present al- 
ready in antiferromagnetically coupled two-layers. [9] As we have seen, phase diagrams for 
multilayer systems with N > 2 become very complex. Hence, systems with few layers are 
probably a better starting point for detailed investigations on magnetization processes. Such 
experiments could be used to assess magnetic parameters and quality of such systems. Gen- 
erally, magnetization processes and checks for their dependence on magnetic and thermal 
pre-history should be made by applying fields in oblique directions and/or under rotating 
fields. Only then, the behaviour of the ( if ^,7/^,) -phase diagrams can be usefully compared 
with detailed theoretical investigations. 
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Further theoretical work may address models with non-equal constants. Also effectively 
ferrimagnetic systems with odd numbers of layers and with different layer thicknesses 
may be interesting. Some experimental data for magnetization processes in such systems 
exist, e.g. for Co/Cr two-layers with different ferromagnetic layers thicknesses, [14] and for 
odd-numbered multilayers. j9|] Interesting reorientational effects should arise also in magnetic 
fields, which are applied perpendicular or under arbitrary angle to the layer plane. Anti- 
ferromagnetically coupled superlattices may also undergo transitions into perpendicularly 
magnetized states for certain thicknesses of the individual ferromagnetic layers, as observed 

n 

for Co/Cr(001). 12] In such cases the stray field must be taken into account already for 
laterally homogeneous states. 

Concluding we state that the magnetic effects and phenomena discussed in this paper can 
be used for detailed investigations on such aspects of nanomagnetism as interlayer-exchange 
interactions, reorientational, and multidomain processes. 
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